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1.  Introduction 


The  systems  carrying  electric  charges  in  the  form  of  plasma  are  widespread  in 
nature.  Various  US  Army-related  applications  deal  with  different  plasmas.  We  do 
not  dwell  on  any  particular  applications  since  we  have  analyzed  theoretical  issues 
equally  relevant  for  any  of  those  practical  applications. 

Basically,  we  address  the  general  problems  related  to  thermodynamics  of  plasma, 
including  equilibrium  and  stability  conditions  for  macroscopic  systems  containing 
components  in  a  plasma  state.  Of  course,  these  general  issues  have  been  addressed 
and  revised  in  tens  of  classical  textbooks,  and  certainly  they  will  be  addressed  and 
revised  in  hundreds  of  textbooks  to  come.  This  keeps  happening  for  2  main 
reasons.  First,  the  fundamentals  of  any  discipline  are  full  of  striking 
inconsistences  and  contradictions.  These  striking  contradictions  will  never  be 
fully  avoided,  although  they  stay  dormant,  sometimes,  for  decades  and  even 
centuries.  Second,  the  number  of  promising  and  useful  applications  keeps 
growing  and  requires  permanent  revisions  of  fundamentals. 

Thus,  we  do  not  address  all  fundamental  problems  and  contradictions;  this  is 
simply  impossible.  Instead,  we  discuss  the  bare  minimum  of  the  thermodynamical 
tools — only  those  tools  that  are  absolutely  unavoidable  when  dealing  with  the 
problems  of  equilibrium  and  stability  of  systems  containing  gaseous  plasmas. 

There  are  2  different  approaches  in  thermodynamics,  each  having  its  own 
advantages  and  disadvantages.  Josiah  Willard  Gibbs  made  the  major  contributions 
to  each  of  them  in  his  seminal  treatises. Current  textbooks  and  monograph 
presentations  of  the  thermodynamics  of  plasma  mostly  follow  the  Gibbs^  basic 
principles  of  statistical  thermodynamics  and  their  subsequent  development. 
Landau  and  Lifshitz^  summarize  those  developments. 

The  approach  of  Gibbs, ^  based  on  the  concept  of  heterogeneous  systems  and 
usage  of  variational  principles,  did  not  get  enough  attention.  Here  we  apply  the 
Gibbs  variational  method  (as  we  understand  it)  for  the  analysis  of  the  simplest 
possible  systems  containing  ionized  gases. 


2.  Gibbs  Method  in  the  Integral  Form 

As  per  the  Gibbs  general  methodology,  based  on  the  concept  of  heterogeneous 
systems,  we  have  to  choose  the  energy  and  entropy  functionals  and  specify  the 
macroscopic  degrees  of  freedom  for  the  system.  Consider  a  closed  vessel 
containing  a  gas  or  liquid,  the  material  particles  of  which  carry  charges. 
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Typically,  the  net  charge  of  the  system  is  close  to  zero.  However,  for  the  sake  of 
brevity,  we  assume  that  all  the  charges  are  the  same;  for  instance,  electrons  or 
negatively  charged  ions. 

We  use  the  Eulerian  approach  to  describe  liquid  continuum  media.  Fortunately, 
the  Eulerian  description  is  the  most  convenient  when  dealing  with  both  the 
thermodynamics  of  liquids  and  gases  and  also  with  electromagnetism. 

2.1  The  Functional  for  One-Component  Plasma 

Let  Q  be  the  volume  inside  the  closed  vessel.  Let  p{z)  be  the  mass  densities  per 
unit  volume  of  the  gas  (liquid)  under  study.  Let  <7  be  the  charge  density  per  unit 
mass.  The  charge  density  q  per  unit  volume,  then,  is  given  by  the  relationship 


q  =  (7p.  (1) 

Let  M  and  Q  be  the  total  mass  and  charge  of  the  liquid,  respectively.  They  are 
expressed  by  the  following  integrals  over  the  domain  occupied  by  the  charged 
gaseous  substance: 


M=\d^p{z),  Q  =  \dOq{z)  =  aM.  (2) 

Q  Q 

Let  the  77  (z)  be  the  entropy  density  per  unit  mass  of  the  electric  liquid,  whereas 
S  is  the  total  entropy  of  the  system;  then  we  get 

S=\dapq{z).  (3) 

Let  the  c(z)  be  the  spatial  distribution  of  the  internal  energy  density  per  unit 
mass  of  the  electric  liquid,  whereas  U  is  the  total  internal  energy  of  the  system, 
as  follows: 


U=\dVlpe{z).  (4) 

Similarly,  we  introduce  be  the  free  energy  density  y{z)  per  unit  mass  of  the 
electric  liquid,  whereas  F  is  the  total  internal  energy  of  the  system  and  other 
thermodynamic  functions. 

The  thermodynamics  of  the  system  will  be  completely  defined  when  the  internal 
energy  density  e  is  given  as  a  function  of  the  mass  density  p  and  the  mass 
entropy  density  q ,  as  follows: 
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e^e[p,rj)  . 


(5) 


Given  the  function  e{^p,T]),  there  is  no  need  to  separately  choose  any  additional 
equations  of  state  (EOSs);  like,  for  instance,  the  pressure  p  or  the  absolute 
temperature  T  as  functions  of  2  macroscopic  thermodynamic  parameters.  This 
information  is  already  contained  in  the  function  ei^p^rf)  and  can  be  extracted 

from  this  function  with  the  help  of  the  formula 

=  =  (6) 

op  Ot] 

but  not  the  other  way  around.  Given  the  function  p{p,T])  we  cannot  recover 
from  it  the  internal  energy  density  function  e{yp,rf).  That  is  why  the  function 
eiyp,'q)  can  be  coined  a  complete  EoS,  whereas  piyp,r])  or  T{^p,r])  are 
incomplete  EoSs.  So,  there  are  several  incomplete  EoSs,  none  of  which  contains 
as  much  information  as  eiyp,rj).  However,  several  incomplete  EoSs  can  contain 
altogether  as  much  information  as  one  complete  EoS.  The  advantage  of  the 
incomplete  EoSs  is  implied  by  the  fact  that  they  can  be  more  easily  extracted  from 
physical  experiments. 

One  can  ask,  “Is  there  only  one  complete  EoS  for  given  liquid  substance?”  The 
answer  is  negative.  In  particular,  the  free-energy  density  function  y/{^p,T^ 
contains  as  much  information  as  the  complete  EoS  function  e{^p,rj).  Thus, 
^  (yC>,r)  is  also  a  complete  EOS  on  its  own  right. 

The  statistical  thermodynamics  provides  a  theoretical  procedure  allowing  us  to 
extract  the  complete  EoS  y/{^p,T)  from  data  related  to  the  spectrum  of  existing 
energy  levels  of  the  system.  Calculation  of  the  EoS  still  relies  on  an  additional 
hypothesis  that  is  not  that  easy  to  justify.  Fortunately,  some  of  the  important 
qualitative  facts  can  be  extracted  without  full  knowledge  of  the  spectra. 

The  macroscopic  Gibbs  approach,  based  on  the  concept  of  heterogeneous 
systems,  does  not  suggest  any  procedure  for  calculating  the  complete  EoS  from 
theoretical  reasoning.  Basically,  with  this  approach  the  EoS  should  be  extracted 
from  specially  designed  experiments  dictated  by  the  basic  thermodynamic 
principles,  physical  intuition,  or  a  combination. 

The  peculiarity  of  the  heterogeneous  systems  with  plasma  components  consists  in 
the  necessity  of  taking  into  account  the  nonlocal  electrostatic  energy  as  an 
essential  addition  to  the  classical  internal  energy. 

The  classical,  local,  and  additive  internal  energy  can  be  presented  in  the  form  of  a 
single  integral  over  the  domain  Q ,  occupied  by  the  substance  under  study.  At  the 
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same  time,  the  electrostatic  energy  is  different;  it  is  nonadditive  and  nonlocal. 
Often  it  is  presented  in  the  form  of  a  double  integral  over  the  domain  occupied  by 
the  charged  substance.  The  total  electrostatic  energy  of  the  system  we 

postulate  in  the  following  form  of  the  integral, 


E.,..  =  ipapn- 


(7) 


with  a  singular  core  (i.e.,  within  the  domain  Q  the  integrand  in  Eq.  7  approaches 
infinity). 

The  singularity  of  the  integrand  creates  a  variety  of  technical  difficulties.  To 
avoid  them  the  electrostatic  energy  can  also  be  presented  in  the  form  of  a  single 

integral  of  the  electrostatic  field  E(^z)  without  any  singular  core. 


OTT  e 

Space 


(8) 


but  the  integral  now  spans  over  the  whole  space,  not  just  part  of  it,  occupied  by 
the  ionized  substance.  The  field  components  E^i^z)  are  the  gradients  of  the 
electrostatic  potential  (l>{z) . 


£,(z)  =  -V,<z).  (9) 

The  theoretical  electostatics  were  built  based  on  the  analogy  with  the  earlier 
developed  theory  of  self-gravitating  celestial  bodies.  Still,  many  models  of 
electrostatic  systems  are  almost  the  same  as  the  models  of  the  systems  with  self¬ 
gravitation.  For  instance,  the  total  gravitational  energy  E^^^^  of  a  self-gravitating 
body  is  typically  postulated  in  the  form 


(10) 


We  see  that  the  electrostatic  and  gravitational  energy  are  very  similar  but  differ  by 
the  sign.  The  difference  in  the  signs  makes  the  physics  of  self-gravitating  and  self- 
repelling  electrostatic  systems  completely  different,  although  the  mathematical 
tools  used  in  both  disciplines  are  almost  identical. 
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In  the  following,  we  neglect  the  gravitational  energy.  We  assume  that  the  total 
energy  of  the  system  comprises  the  total  internal  energy  internal  and 
total  electrostatic  energy  given  by  the  integrals 


Kt  =  Ki+Kiec  =  \d^pe{p,il)  +  ^\d^\da 


Q 


^p{z)p[z) 


(11) 


Q 


Q 


2.2  Basic  Gibbs  Variational  Principle 

Following  Gibbs,  we  postulate  that  the  thermodynamic  and  mechanical 
parameters  of  the  equilibrium  configuration  deliver  the  minimum  (more  precisely, 
stationary  value)  to  the  total  energy  E,^,,  given  by  Eq.  11.  When  talking  about 
minimum,  we  must  specify  the  mechanical  and  thermodynamical  degrees  of 
freedom  of  the  system  under  consideration.  Quite  often,  it  is  done  by  specifying 
allowable  infinitesimal  variations  of  certain  parameters.  We  postulated  that  for 
our  system  the  virtual  variation  of  the  mass  density  5p  and  entropy  St]  must 
keep  the  following  integral  amounts  fixed: 

\dCLp{z)  =  M,  (12) 

n 


and 


jrfQ77(z)  =  5.  (13) 

n 

Using  the  standard  method  of  the  Lagrange  indefinite  multipliers,  we  arrive  at  the 
following  unconditional  variational  problem  for  the  functional  O  : 

O  =  \dap \e{p,r^)-h-@^^']  +  —  f^/Q fi/Q*  ^  ,  (14) 

n  ^  Cl  Cl 

where  A  and  0  are  the  indefinite  Lagrange  multipliers. 

2.3  First  Variation  and  the  Equilibrium  Equations 

The  first  variation  of  the  functional  O  reads 
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^  =  jdQ 


/  L 


Q 


p(er^-%)5ri 

Indeed,  the  functional  <I)  can  be  presented  as  sum  of  2  terms,  as  follows: 


(15) 


where 


I^=\dVlp[e{p,r])-®ri-A\ 


(16) 


(17) 


and 


pU)p(z) 


n  n 


z  -z 


Varying  the  expressions  Eqs.  17  and  18,  we  get,  respectively. 


(18) 


SIn=jdQ  i^{pe)^-®ri-KYp  +  p[e^-®)Sri^  (19) 


and 


dI^^cj^\dCLdp{z)\d^* 

Q  Q 


p(z) 

I- 


(20) 


The  denominator  2  in  the  relationship  Eq.  20  disappeared  because  p  in  the 
integrand  of  Eq.  18  appears  twice  as/7(z)  andp(^z*) . 


Adding  the  relationships  Eqs.  19  and  20  we  arrive  at  the  required  relationship 
Eq.  15  for  the  first  variation  of  ^ . 

Separating  the  independent  variations  5p  and  St]  in  Eq.  15,  we  arrive  at  the 
following  conditions  of  the  bulk  equilibrium: 


and 


(yC>e)^  -  ©rj  +  a^^dCl* 

n 


p{z) 


=  A 


(21) 


Approved  for  public  release;  distribution  is  unlimited. 

6 


(22) 


e^=@. 

Thus,  we  arrive  at  the  system  of  4  equations  describing  equilibrium 
configurations:  3  integral  Eqs.  12,  13,  and  21,  and  1  algebraic  Eq.  22,  with  respect 
to  4  unknowns:  2  spatial  function  yC»(z),?7(z),  and  2  unknown  constants  A  and 
0.  No  boundary  conditions  are  necessary  for  this  system. 

Let  us  introduce  an  electrostatic  potential  (l>{z),  defined  as 

^(z)  =  crfjQ*  .  (23) 

n  I 

We  then  can  rewrite  the  relationship  Eq.  21  as 

e{p,ri)-®ri  +  pep{p,T^)  +  (7(l)  =  h.  (24) 

Using  thermodynamic  identities,  Eq.  6,  we  can  rewrite  the  relationship  Eq.  24  as 

e(^p,ri)-Tri  +  —  +  (7<l)  =  K.  (25) 

P 

The  relationship  Eq.  25  can  also  be  rewritten  in  terms  of  the  free  energy  density 
xp  =  e-rjT  as 


y/  +  —  +  <7^  =  A .  (26) 

P 

The  relationships  Eqs.  25  and  26  are  equivalent  from  the  standpoint  of 
mathematics;  however,  depending  on  circumstances,  one  or  another  form  appears 
to  be  more  technically  convenient. 

In  the  absence  of  electric  field  the  relationship  Eq.  25  implies 

e-Tr]  +  ^  =  A.  (27) 

P 

The  combination  on  the  left-hand  side  of  Eq.  27  is  the  specific  value  (per  unit 
mass)  of  the  so-called  Gibbs  thermodynamic  potential,  or  the  Grand 
thermodynamic  potential  of  one-component  liquid  substance.  For  heterogeneous 
systems  containing  2  one-component  phases,  this  quantity  also  plays  the  role  of 
the  chemical  potential  p  of  the  phases.  In  other  words,  the  full  equilibrium  in  the 
heterogeneous  system,  containing  2  one-component  phases — not  only  the 
pressures  and  temperatures  of  the  phases — should  be  equal  ( Pi  =  P2  ’  ^1  =  ^2 ) »  but 
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the  chemical  potentials  of  the  phases  (//j  =  iJ.j)  should  also  be  equal.  We  do  not 
proceed  with  this  discussion  in  this  paper;  it  will  be  done  later.  The  interested 
reader  should  refer  to  the  Grinfeld  monograph.^ 

The  relationship  Eq.  25  shows  that  when  dealing  with  electrostatic  forces  the 
additional  term  cr^,  responsible  for  electrostatic  interaction,  should  be  included  in 
the  equilibrium  equation  and  in  the  chemical  potential  //.  This  fact  was  known  to 
Gibbs,  ^  who  established  it  from  a  different,  less  formal,  and  more  intuitive 
reasoning.  The  advantage  of  the  reasoning,  presented  here,  is  that  it  allows  us  to 
proceed  with  calculation  of  the  second  variation  and,  thus,  to  get  the  toll  for 
analysis  of  stability  conditions  and  not  just  the  conditions  of  equilibrium.  This  is 
something  that  Gibbs  was  not  able  to  accomplish  in  his  time. 


3.  Second  Variation  and  Stability  Conditions 

The  main  instrument  in  investigation  of  stability  of  heterogeneous  systems  is  the 
second  energy  variation.  If  the  second  variation  is  negative  for  some  allowable 
variations  of  the  thermodynamics  degrees  of  freedom,  we  can  conclude  that  the 
equilibrium  configuration  under  study  is  unstable. 

By  varying  the  relationship  Eq.  15  one  more  time  in  the  vicinity  of  equilibrium 
configuration,  we  arrive  at  the  following  formula  of  the  second  variation: 


d^^{a,h)  =  \d^(Ay  +  lA^^ah  +  + 

Q 

P  P  a(z)a(z*) 

Q  Q,  F  ~  ^ 


(28) 


where 


a{z)  =  Sp,  h{z)  =  Sti, 

(29) 

^PP  ~  ’  ^pr?  ~  ~  ^^'7'? ' 

For  stability,  the  integral  quadratic  form  5^<b{^a,h)  should  be  nonnegative  for  the 
arbitrary  variations  satisfying  the  linear  bulk  constraints. 

ji/Qa(z)  =  0,  ji/Q/i(z)  =  0.  (30) 


Approved  for  public  release;  distribution  is  unlimited. 

8 


3.1  Spectral  Analysis  of  the  Second  Variation 


Consider  the  minimum  of  the  second  variation  Eq.  28  under  the  isoperimetric 
constraints  of  Eq.  30  and  the  normalization  condition 

\dO£i\z)  =  \.  (31) 

n 

As  before,  we  arrive  at  the  unconditional  minimization  of  the  functional 

+ 

f  oi  1  2  • 

J  dO.  — V-  -cca-  Bn  -  Aa 

a  I?-?  I 

In  the  relationship  Eq.  32,  a,  P,  and  X  are  the  indefinite  multipliers  associated 
with  the  constraints  Eqs.  30  and  31. 

Eor  the  first  variation  of  this  functional  we  get 


=  Jr/Q 

n 


-SJ{^a,h)  = 


V  '  Q  r  ^  / 

Thus,  we  arrive  at  the  nonuniform  linear  bulk  equations 


(^pp  +  +  jrfQ 


1 

=  —a 

I  I  O 

z  -z  ^ 


.(33) 


(34) 


and 


+  =  (35) 

We  arrived  at  5  equations,  2  Eqs.  30,  plus  Eqs.  31,  34,  and  35,  and  5  unknowns, 
a,h,a,p,X.  Consider  a  solution  of  this  system  to  be 


a  =  a°  ,h  =  h°  ,a  =  a°  ,p  =  P°  ,sjLidX  =  X° ,  (36) 

where  the  superscript  “  o  ”  refers  to  the  equilibrium  values  of  the  corresponding 
parameters. 

So,  by  definition,  we  get  the  following  relationships: 
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(37) 


\dQa\z)  =  \, 

J^iQ/i°(z)  =  0,JjQfl°(z)  =  0,  (38) 

Q  Q 

{A„-^y  +  A„h’  +  ida-  =  ia-  ,  (39) 

n  2 


and 


+  ■  (^0) 

Let  us  consider  Eqs.  38,  39,  and  40  as  a  system  of  4  linear  equations  with  4 
unknowns  a,h,a,and/3.  This  system  is  not  only  near  but  also  uniform. 
Therefore,  it  always  has  a  trivial  solution  a  =  h  =  a  =  f3  =  0.  This  solution  is  of  no 
interest  to  us  since  in  view  of  Eq.  37  we  need  a  solution  oft  a°  ^  0.  So,  we  have  to 
have  a  situation  in  which  the  system  of  Eqs.  17a,  18,  19  has  multiple  solutions. 
This  is  possible  only  for  some  special  values  of  y.  We  call  them  spectral  values. 
In  the  following  it  is  evident  that  the  spectral  values  play  a  special  role  in  the 
analysis  of  stability  of  equilibrium 

For  a  spectral  value  X° ,  let  us  calculate  the  corresponding  no-vanishing  values  of 
the  unknowns  a  ,h°  ,a° ,  and /3° .  Let  us  multiply  Eq.  39  by  a  ,  Eq.  40  by  h° ,  and 
integrate  the  resulting  relationships  over  the  volume  Q.  Then,  using  Eq.  17,  we 
get  the  relationships 


(z)  +  {z)  +  jdQjdQ* 

Q  Q  Q  Q 


\z-z  \ 


(41) 


and 


+  jdQA^^h°^  =  0 .  (42) 

n  n 

Adding  the  relationships  Eqs.  41  and  42,  we  get  the  relationship 


jdn(A^^a^  +2A^^h°a  +  +  jdnjdQ^ 

n  n  n 


a°(z)a°( 

y) 

-i- 

z  -z 

=  r  .  (43) 
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Comparing  Eq.  28  with  Eq.  43,  we  get 


5^^[a,K^)  =  r.  (44) 

In  words,  the  relationship  Eq.  44  says  that  the  spectral  values  A  and  the 
associated  nontrivial  solutions  of  the  system  Eqs.  38,  39,  and  40  are  equal  to  the 
extrema  of  the  second  energy  variation  of  the  system.  Thus,  we  get  the  following 
necessary  conditions  of  stability.  For  stability  all  the  solutions  of  the  system 
Eqs.  29,  31,  34,  and  35,  all  the  values  of  X°  must  be  nonnegative. 

We  can  reformulate  this  statement  as  the  following: 

For  stability,  all  the  spectral  values  X  of  the  linear  uniform  system 

\dQa  (z)  =  0,  \dah°  (z)  =  0,  (45) 

Q.  Q 


{A^,-A)a‘+AJ,‘+yn- 

Q. 


(46) 


and 


+  (47) 

with  respect  to  the  unknowns  a,h,a,  and  must  be  nonnegative.  In  the 
following  we  will  be  calling  this  statement  the  Stability  Principle  of  Nonnegative 
Spectrum. 

3.2  Thermodynamic  Inequalities 

The  problem  of  stability  and  the  stability  principle  of  nonnegative  spectrum  has 
an  immediate  relation  to  the  problem  of  thermodynamic  inequalities.  For  the 
analysis  of  thermodynamical  inequalities  it  is  convenient  to  rewrite  Eq.  46  as  the 
following  pair  of  equations: 


(A„-Ay+A^,h’+f~a’=0  (48) 

VV.f  +47ra°  =  0.  (49) 

Consider  a  shortwave  spectrum  of  the  system  Eqs.  47^9  for  which 
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a°  ^(5°  ^  0,  a°  =  ,  f  =  . 


(50) 


Then,  Eqs.  47^9  imply  the  following  system: 


(51) 

a^-4^++,h-^-'=o. 

(52) 

and 


J  +  AttA  =  0 


(53) 


By  eliminating  the  scalars  H  and  J  in  Eqs.  51-53,  we  arrive  at  the  single  linear 
uniform  equation 


-2  A^  ^ 

k 

pn 

AjTj  J 

A  =  0 


with  respect  to  the  remaining  constant  A. 

Equation  54  always  has  the  following  trivial  solution: 


(54) 


A  =  0.  (55) 

The  nontrivial  solutions  exist  only  when  the  expression  on  the  brackets  in  Eq.  54 
vanishes.  This  obvious  fact  leads  us  to  the  following  formula  for  the  spectral 
values  X  : 


X  = 


rji] 


■  +  471 


-2 


Thus,  we  arrive  at  the  inequality 


pp 


el 


+  4k 


>0. 


If,  instead  of  the  normalization  condition  Eq.  31,  we  use 


\dQh\z)  =  \, 


then  instead  of  the  system  Eqs.  51-53,  we  get  the  system 
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(56) 


(57) 


(58) 


AA  +  AH  +  J=Q. 

PP  Pri 


and 


k\  J  +  A7rA  =  0. 


System  Eqs.  59-61  leads  us  to  somewhat  different  spectral  values  of  X. 


A,: 


-2 


-2 


Thus,  we  arrive  at  the  thermodynamic  inequality 


-2 


X^nn  +  4^ 

PP 


-2 


>0. 


At  last,  let  us  consider  the  case  of  the  normalization  condition 

|jQ(a"  +  /eV)  =  l. 

n 

In  this  case,  we  arrive  at  the  spectral  system 


K-2)a+A„«  +  2=0, 


and 


J  +  4;7rA  =  0 , 


which  can  be  rewritten  in  the  following  matrix  form: 


(59) 

(60) 

(61) 

(62) 


(63) 

(64) 

(65) 

(66) 

(67) 
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(68) 


pt! 

PP 

-2 


At 


A7r\k 

00,  the  Eq.  68  reads 


pt] 

0 


0 

1 

-1 


r  A' 

'o' 

H 

0 

[y 

0 

^pt] 

o' 

'A' 

‘O' 

\p-^ 

^pt] 

1 

H 

= 

0 

0 

0 

-1 

J 

0 

(69) 


leading  to  the  spectrum  (secular)  equation 


/l"- 


\t]  ^  ^PP_  ^  _|_  \v^PP  ^PV 


2  '*7'/  '  ‘^PP 

K 


=  0. 


(70) 


Equation  71  has  the  following  discriminant  A  : 


=(4,  +  -4;.=  (\,A„-Al,)  = 

+Ak-a!-^ 


(71) 


Obviously,  the  discriminant  A  is  always  positive,  and  therefore  the  spectrum 
values  A.  are  real.  They  will  be  nonnegative  provided  the  following  relationships 
are  satisfied: 


and 


A„+k^A>Q. 

rji!  PP 


(72) 


(73) 


Since  k  in  the  inequality  Eq.  73  can  be  an  arbitrary  real  constant,  the  positiveness 
of  the  Eigenvalues  A  implies  2  independent  inequalities: 


(74) 


Thus,  in  the  asymptotics  k  we  arrive  at  the  3  classical  thermodynamic 
inequalities  of  Eqs.  72,  73,  and  74. 
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4.  Conclusion 


We  established  general  equations  allowing  analysis  of  the  equilibrium 
configurations  of  the  systems  containing  electrically  charged  constituents.  For  the 
sake  of  simplicity,  technical  transparency,  and  brevity,  we  limited  ourselves  to  the 
systems  containing  charges  of  a  one  sign.  Our  approach  was  based  on  the 
variational  principles  of  Gibbs,  which  are,  in  turn,  based  on  the  concept  of 
heterogeneous  systems.  The  deduction  of  the  equation  of  equilibrium  is  based  on 
the  calculation  of  the  first  energy  variation  of  the  functionals  with  isoperimetric 
constraints.  We  arrived  at  the  system  of  4  equations  describing  equilibrium 
configurations:  3  integral  Eqs.  12,  13,  and  21  and  one  algebraic  Eq.  22  with 
respect  to  4  unknown:  2  spatial  functions  p{z)  and  2  unknown 

constants,  A  and  0.  No  boundary  conditions  are  necessary  for  this  system. 

We  then  established  the  necessary  conditions  of  thermodynamic  stability  of  the 
corresponding  equilibrium  configuration.  Our  approach  is  based  on  the  derivation 
and  analysis  of  the  second  variation,  which  is  given  by  the  relationship  Eq.  28, 
which  uses  the  notation  Eq.  29.  The  second  variation  appears  in  the  quadratic 
integral  form,  which  should  be  analyzed  in  conjunction  with  2  linear  integral 
constraints  of  the  isoperimetric  type  in  Eq.  30. 

We  demonstrated  how  the  concept  of  stability  can  be  applied  to  the  classical 
problem  of  thermodynamic  inequalities.  We  also  established  the  novel 
thermodynamic  inequalities  Eqs.  57  and  63,  which  generalize  the  classical 
thermodynamic  inequalities  of  Gibbs  for  the  charges  liquids  or  gases. 
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